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Abstract. The high precision gravity measurements to be made by recently launched (and recently 
approved) satellites place new demands on models of Earth, atmospheric, and oceanic tides. The 
latter is the most problematic. T1 e ocean tides induce variations in the Earth’s geoid by amounts 
that far exceed the new satellite sensitivities, and tidal models must be used to correct for this. Two 
methods are used here to determine the standard errors in current ocean tide models. At long wave- 
lengths these errors exceed the sensitivity of the GRACE mission. Tidal errors will not prevent the 
new satellite missions from improving our knowledge of the geopotential by orders of magnitude, but 
the errors may well contaminate GRACE estimates of temporal variations in gravity. Solar tides are 
especially problematic because of their long alias periods. The satellite data may be used to improve 
tidal models once a sufficiently long time series is obtained. Improvements in the long-wavelength 
components of lunar tides are especially promising. 


1. Introduction 

With the recent launches of the CHAMP and GRACE satellite missions and the 
approval of the GOCE mission, we find ourselves entering a new era in satellite 
gravimetry. The high precisions and increased sensitivities of these missions place 
new demands on methodologies and models. This paper examines the role of tides 
and tide model errors, with special emphasis given to the GRACE low-low satellite- 
to-satellite tracking mission (Tapley and Reigber, 2000). 

Tides have been a key component of space geodesy since the mid- 1960’s, but 
their role in satellite gravimetry has been confined primarily (although not exclu- 
sively) to determining long-period tidal perturbations in satellite orbits or to esti- 
mating certain tidal parameters from such perturbations. Determining long-period 
orbit perturbations will conlinue to be of utmost importance, in these new missions 
as well as in other high-precision applications (e.g., Pavlis and Iorio, 2002). But 
in this new era, attention will also turn to two new aspects (cf. Schrama, 1996). 
(1) The continuous GPS tracking of satellites now allows short-period tidally in- 
duced orbit perturbations to be studied, necessitating models with complete sets 
of spherical harmonic coefficients, not limited to a few select orders. (2) Some 
of the new data types allow a direct instantaneous measurement of the gravity 
effects of tides (embedded, of course, in a background of gravity signals from other 
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anomaly sources). Unlike long-period orbit perturbations, these direct effects can 
be localized near the causative masses. Furthermore, the contamination of gravity 
field measurements (or estimates) from inadequately modeled tidal masses will be 
an important concern. In this paper we consider the gravitational effects of tides 
and examine how errors in tidal models may affect the new satellite gravimetry 
missions. Space precludes a complete discussion, so we focus on oceanic tides, as 
opposed to Earth tides or atmospheric tides, because errors in the former are much 
larger than errors in the latter. 


2. Ocean Tides and Gravity 

To understand the extent to which new gravity missions are directly affected by 
ocean tides, we use a global model to evaluate ( 1 ) the tidally induced perturbations 
in gravity (or the geoid) and (2) the associated signals that can be expected in 
satellite measurements. Later sections consider corresponding errors in present-day 
models. 

2.1. Gravity Perturbations Induced by Ocean Tides 

Suppose the ocean tidal elevations are expanded in normalized complex spherical 
harmonics y„ m (0,(p) as 

«e,9,0 = £W0C(e.<p). d) 

n } m 

The coefficients z nm {t) vary with tidal periodicity. At satellite altitudes the gravita- 
tional potential of this tide is given by (e.g., Lambeck, 1988) 

= 4nGop w £ (-)" + ' z_(i) C(8,<P) (2) 

rt, r 

where p w is the mean density of seawater (approximately 1035 kgm -3 ), a the 
radius of the Earth, G is the Newtonian constant, and k' n are loading Love numbers. 

The corresponding induced variations in the geoid are given by (e.g., Schrama, 
1996; Wahr, 1998) 

8/V(0,cp,r) = 3(p*/p e ) V Znm{t) C(0><P)> (3) 

where p e is the mean density of the Earth. Figure 1 shows the amplitudes of both ^ 
and 57V for the principal semidiurnal lunar tide M 2 . 

* Here ‘direct’ is used in contrast with the indirect tidal gravity effect manifested by orbit pertur- 
bations. Schrama (1996) prefers ‘indirect’ for both and reserves ‘direct for the primary astronomical 
tidal forces of the sun or moon. The latter are exceedingly well known and not discussed in this paper. 
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Figure I. Amplitude of the M 2 < >cean tide (top) in cm and the corresponding induced geoid per- 
turbation (bottom) in mm. The dean model is based primarily on Topex/Poseidon data. The geoid 
perturbation follows from Eq. (3) 


For reference the tidal perturbation in the dimensionless Stokes coefficients C„ 
are related to the elevation c oefficients by 

e r ( f \ _ ^Pi v (l4-& n ) _ , , 

oCnm\t) — (n . Znm{t )• ' 

ap e (2n+ 1 ) 

At radius r = a the tidal perturbations in the gravity field —3 U /3 r are 


8g(0,(p,/) =47tGp w X( l + k 'n) 


2 n+ 1 


Znm ( f )y;(e,(p). 


For sufficiently large n the gravity perturbation is related to the tidal height simply 
by 


lim 8 g nm = 2nGp w z nm « (0.0434 mGal/m) z nm , 
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which leads to the Bouguer plate approximation (8 g ~ 0.0434Q used by Knudsen 
and Andersen (2002). The approximation works well except for very small n. 

As Figure 1 shows, the tides perturb the geoid by amounts well below 5 cm. 
The largest perturbation (for M 2 ) is 3.7 cm in the equatorial eastern Pacific. Geoid 
signals of this magnitude are about an order of magnitude smaller than present 
uncertainties in the geoid over the ocean (Lemoine et al., 1998), but they are 
well within the anticipated precisions of GRACE and GOCE (for GOCE, see ESA 
[1999] Table 8.4). 

2.2. Ocean Tide Signai s in GRACE and GOCE 

The gravitational signals of ocean tides relevant to the GRACE and GOCE mea- 
surements can be predicted with a global ocean tide model. Cheng (2002) recently 
employed linear perturbation methods to compute the anticipated spectrum of tidal 
signals in GRACE data over a frequency range of 0-150 cycles/day. We here 
show an example in the time domain based on a numerical integration of GRACE 
ephemerides (including initial state adjustments) at altitude approximately 500 
km. The relevant forcing from ocean tides is computed using the the four largest 
constituents of the GOT00.2 ocean model (an update to Ray [1999]) complete to 
spherical harmonic degree 60. 

Figure 2a shows the expected ocean tide signal in GRACE’S range-rate mea- 
surement for an arbitrary 1 20-minute arc segment, or about 1 .5 revolutions. Most 
of this signal should be observable with GRACE’S l-/im/s precision, especially 
given that these signals are temporally coherent. 

As an alternative to employing numerical integration or linear perturbation meth- 
ods, it is sometimes useful to approximate these range rate signals by employing 
the following simplistic method. For a perfectly circular and polar orbit, the along- 
track acceleration is given by 0- V£7 and the velocity thus by ( \ j r) j (dU /o)0) dt . 
Within the integrand we approximate 0 as simply (GM/r 3 ) 1 / 2 , with M the Earth’s 
mass, so that the along-track velocity integrates to v(r) r; (r/GM)^ 2 U{t). The 
range rate p is then taken as the difference v(f) — v(/ + A/), with At the time lag 
between the two satellites, <>r 

p(0«(GM/r 3 ) 1 / 2 (G(n-(/(t + Ar). (6) 

Of course, this simple approximation fails to account for the dynamically resonant 
terms in range rate at 0, 1, and 2 cycles/rev; these terms, which in fact are dominant 
for GRACE (Cheng, 2002) can only be recovered by integration of the equations 
of motion. Nonetheless, the method usefully captures the shorter period signals in 
range-rate, as can be seen in Figure 2b. The difference between Figure 2a and 2b 
is predominantly a sum of two sinusoids of frequency 1 and 2 cycles/rev. 

The planned GOCE satellite will measure directly the nine components of the 
second-order gravity gradient tensor. Figure 3 shows the predicted tidal signal in 
the second-order radial gradient, which is generally the largest component (Schrama, 
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Figure 2. Predicted ocean tide signals in the GRACE range rate measurements for an arbitrary 120 
minutes of a 1-day arc, computed (a) by numerical integration with a complete dynamical force 
model and (b) by the approximation in Eq. (6). The difference between these curves is caused 
primarily by 1- and 2-cycle/rev resonant terms, not accounted for in (6). 


1996). The arc is the same 120-minute arc used above, but artificially lowered to 
250-km altitude. The signal is computed simply as the gradient of the potential U in 
(2). As should be expected, a gradiometer is sensitive to relatively short wavelength 
components of the tidal mass field. In fact, the signal in Figure 3 tends somewhat to 
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Figure 3. Predicted ocean tide signals in the GOCE gradiometer measurements of d 2 U /dz 2 where U 
is given by Eq. (2), taken to degree and order 150 for the four largest tidal constituents. The arc is 
identical to that used for Figure 2 except artificially lowered to 250 km altitude. Vertical axis units 
are /xE, or 10” 15 s" 2 . 
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mimic the tidal elevation field directly below the satellite. For example, the sharp 
jump near t = 25 m occurs as the satellite path passes from land to a coastal region 
of large tide. Unlike the case for GRACE, the tidal signals in GOCE are more 
than an order of magnitude smaller than the instrument measurement precision, 
which for the gradiometer is a few mE/Hz 1 / 2 (where 1 Eotvos = 10 9 s 2 ). Within 
the gradiometer measurement band of 0.005-0.1 Hz (ESA, 1999), tidal signals are 
likely to be detectable only by exploiting their temporal coherence. GPS tracking of 
GOCE strengthens the low-frequency measurement precision substantially, which 
further strengthens the tidal sensitivity of the system (see error spectrum below). 

We turn now to a discussion of how well these tidal signals can be modeled as 
an a priori correction to the satellite measurements. 


3. Errors in Present-Day Tide Models 

Comparisons of models with in situ tide gauge estimates suggest that open ocean 
tides are now known to better than 2 cm rms (e.g., Shum et al., 1997; Ray, 1999). 
Errors in coastal regions, where tides have large amplitudes and high wavenumbers, 
are probably an order of magnitude worse. Errors are also substantially larger in 
polar regions above the ±66° latitude limits of Topex/Poseidon data. We require 
here a more definite error model, and we consider two types. First, we simply 
difference some of the better global tide models now available. Since these mod- 
els were all developed from Topex/Poseidon measurements, the resulting error 
model may be somewhat optimistic. Secondly, we adopt the global inverse error 
estimates described by Egbert and Erofeeva (2002). These are based on a Monte 
Carlo method that is somewhat dependent on assumed prior error covariances for 
the altimeter data and for the dynamical tidal equations to which the altimeter data 
are fit. Employing both types of error models is therefore a useful test of consis- 
tency. We concentrate on the M 2 constituent; errors in other tidal constituents are 
generally smaller, with the possible exception of S 2 errors in low latitudes where 
contamination from the S 2 atmospheric tide is apparent in some tide models at the 
level of 1 cm. The error models are described in more detail in Ray et al. (2001). 

Figure 4 shows the geoid amplitude spectrum of the M 2 tide and our two dif- 
ferent error estimates. These spectra are computed from the complex Stokes co- 
efficients 5 C nm (t) (or from their estimated errors) following Eq. (4). The degree 
amplitude spectrum in the geoid is 

{ n \ 1/2 

SN n = al X |SCm>| 2 I > (7) 

\m=-« ) 

where the overbar denotes the rms over a tidal cycle (12.4 h for M 2 ) in each com- 
plex component of 5 C nm (t) This form of 6N n is consistent with Eq. (4) of Ray et 
al. (2001) and is also consistent with the geoid spectra of Wahr et al. (1998). 
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Figure 4. Degree amplitude geoitl spectrum of the M 2 tide (thin solid line) and two estimates of its 
error (thin broken lines, the inverse model being the lower curve at low degree). These curves repre- 
sent rms spectra over a single tidal cycle (12.4 h). Heavy solid lines denote anticipated sensitivities 
of three multi-year space geodetic missions, as estimated by Visser et al. (2002). Tidal constituents 
other than M 2 would fall below these M 2 curves. 


In Figure 4 the M 2 error curves are the broken lines, one being the Monte Carlo 
inverse model, the other being the mean of the differences among models TPXO.5, 
GOTOO.2, and NA099 (the latter from Matsumoto et al., 2000). The inverse error 
is the lower curve at small degree n and the higher at large n. The two error curves 
are reasonably consistent and generally within a factor of 2. The bold lines in 
Figure 4 represent the sensitivities of three different space geodetic missions after 
multiple years of data. As expected, the M 2 signal is well above the sensitivities of 
all three satellite missions, which suggests that high-degree tidal models, perhaps 
to at least degree 70, are required for GRACE and GOCE. Figure 4 also shows 
that present-day M 2 errors exceed the GRACE sensitivity below degree 50 or so. 
The errors also exceed the CHAMP sensitivity below degree 8 or 10. They do not, 
however, exceed the expected GOCE sensitivity anywhere. This, and the fact that 
tidal errors are even smaller than the already very small gradient signals shown 
in Figure 3, suggests that ude model errors are unlikely to be a major concern 
for the GOCE mission. It is possible that some large tidal errors in a few coastal 
regions might degrade GOCE analyses, but such errors are far more problematic 
for satellite altimetry than for satellite gradiometry. 

While Figure 4 suggests that tide models must be improved for (or by) the 
GRACE and CHAMP missions, the temporal mismatch between the tidal and the 
satellite curves (12.4 h vs. several years) complicates matters. The manifestation of 
tidal error in GRACE, for example, depends on issues of sampling and on the de- 
sired estimation timescales (e.g., instantaneous gravity anomalies, monthly means, 
or multi-year means). 
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4. Tidal Errors and GRACE 

The effects of tide model errors on GRACE’S intersatellite measurements are suf- 
ficiently complex that they are difficult to predict without extensive simulations. 
We show one example below. Nonetheless, some insight can likely be gained from 
simple considerations of tidal aliasing. GRACE signals such as those in Figure 
2b will be aliased in ways quite similar to the aliasing in satellite altimetry. One 
must bear in mind, however, that the spectrum of tide-induced signals in GRACE 
is richer than Figure 2b suggests (e.g., Cheng, 2002). GRACE is no altimeter, and 
some of its tidal errors are entangled with orbit errors with potentially complicated 
aliasing. 

4.1. Tidal Aliasing 

The tidal phases observed by an orbiting spacecraft are determined by the angle 
between the satellite orbit plane and the tide-raising body. The orbit plane of any 
high-inclination spacecraft like GRACE precesses so slowly with respect to the 
sun that solar tides are necessarily aliased into relatively long periods — e.g., semi- 
annual, annual, or longer. Such periods are of prime importance to climatological 
studies, making aliased errors at these periods very problematic. (The precession 
rate with respect to the moon is rapid, and aliasing of lunar tides is rarely prob- 
lematic unless one happens to be unlucky.) Because the solar tides that GRACE 
observes will be aliased to long periods, we may predict that errors will be prevalent 
in the GRACE measurements at these same periods, no matter the exact manner in 
which the errors are manifested through dynamical complications in the orbit. 

From the satellite precession rate we may determine the tidal alias periods. The 
alias period of the Si tide is identical to the satellite nodal period with respect to 
the mean sun; the S 2 period is half that. The alias period of Ki is identical to the 
satellite nodal period with respect to the mean equinox; K 2 is half that. The Pi alias 
is identical to the nodal period with respect to a fictitious body moving at twice the 
sun’s mean motion. As of this writing (a few months after launch) the GRACE 
mean Keplerian elements are (a,e,l) = (6860.4km, 0.00256, 89.02°), implying a 
node rate of approximately —0. 132° /day, or —1.1 17°/day with respect to the mean 
sun. The corresponding alias periods for the primary solar tides are: 

K 2 3.74 years Ki 7.48 years 

S 2 161 days Si 322 days 

Pi 171 days 

Some tidal contamination should therefore be expected at the important semi- 
annual period and at the interannual periods of 3.7 and 7.5 years. The Si alias 
is also somewhat problematic, since both Si (being primarily atmospheric) and the 
nearby annual cycle (Sa) have strong temporal variability, causing potential overlap 
of their rather broad spectral peaks. 
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Figure 5. Results of GRACE simulation showing errors in an estimated monthly geoid (in mm) 
caused by tide model errors. Geoid is evaluated to degree and order 22; higher degree expansions 
show larger tidal contamination. Note errors are not confined to the oceans, primarily because of the 
long-period dynamical perturbations to GRACE'S range rate. 
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Figure 6. As in Figure 5 except using the simplistic (non-dynamic) approach of Eq. (6). Errors are 
more localized to the ocean. 


4.2. Simulation of Monthly Gravity Errors 

Because the GRACE project intends to generate monthly mean estimates of the 
geopotential (Tapley and Reigber, 2000), it is of interest to anticipate possible con- 
tamination of this product caused by errors in the tide model adopted for processing 
the GRACE measurements. We have estimated this effect by simulating one (ar- 
bitrarily chosen) month of GRACE range-rate measurements, using the methods 
described by Rowlands et id. (2002). The simulation consists of an inversion for 
a degree- 120 static geopotential in the presence of realistic errors in intersatellite 
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range measurements and in GPS tracking, with and without tide model errors. The 
latter consists of the difference between tide models GOTOO.2 and TPXO.5, to 
degree 60, for constituents M 2 , S 2 , N 2 , K 2 , Oi and K[. The resulting tidal contam- 
ination in the monthly mean geoid estimate is shown in Figure 5. Because errors in 
monthly geoid estimates rise rapidly with degree, the geoid in Figure 5 is evaluated 
only to degree 22. Summing to degree 50 results in fairly similar characteristics, 
although with more orbit stripping and roughly double amplitudes. 

The simulation represented by Figure 5 is computationally expensive, and it 
is therefore of interest to see how well the simplistic (and inexpensive) formula 
(6) can represent the same e rror field (akin to comparing Figures 2a and 2b). This 
involves no ephemeris computations and only a relatively small normal matrix. The 
result is shown in Figure 6. Because this method ignores the dynamical long-period 
effects, there is no track-type stripping and errors are confined primarily to ocean 
regions. The errors in both figures are of comparable magnitude and they align in 
a few locations (e.g., Amazon mouth, Weddell Sea), but the method underlying 
Figure 6 is evidently too simplistic to map the error in detail. 

Several points about Figure 5 are worth noting. The presence of geoid errors 
over both land and ocean is caused by the rich spectrum that tidal forcing induces 
in the range-rate measurements. Localized errors in ocean tide models can induce 
global errors in gravity estimation. Secondly, the amplitude of these geoid errors 
is significant for some anticipated GRACE applications. In particular, any inter- 
pretation of the estimated gravity in terms of oceanic mass motions is considerably 
confounded. Wahr et al. (1998; their Plate 1) estimate annual variations in the geoid 
of about 1 mm and less, based on a general circulation model, and such signals are 
comparable to the errors in Figure 5. Thirdly, the errors in Figure 5 are one monthly 
manifestation and we have no clue to their temporal variations without extending 
the simulations over a considerably longer timespan. From our discussion of alias- 
ing, however, we anticipate that some of these errors have periods of semi-annual 
and longer. The S 2 and K 2 tides are likely to be most problematic for monthly mean 
gravity estimates (Knudsen and Andersen, 2002). 


5. Summary Discussion 

The ocean tides induce variations in the geoid of several cm, well within the an- 
ticipated precisions of the latest space geodetic missions. Yet it is important to 
realize that the geoid errors in present-day gravity models are a few decimeters 
over the ocean (Lemoine et al., 1998), whereas errors in the prediction of tidal 
geoid variations are orders of magnitude smaller. We conclude that missions like 
GRACE and GOCE can improve our knowledge of the geopotential by an order 
of magnitude, or more, even with present ocean tide models. Tide model errors 
become crucial with the extreme precisions required for GRACE to detect temporal 
gravity changes of the sort discussed by Wahr et al. (1998). 
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To the extent that Figure 4 is realistic, GRACE itself offers opportunities to 
improve the long wavelength components of ocean tide models, specifically co- 
efficients at degrees 50 and below, or wavelengths of about 1000 km and longer. 
CHAMP may also help improve a few low-degree coefficients. These improve- 
ments are likely to occur only after multi-year time series have been obtained, 
because many tides are aliased to nearby frequencies. The solar tides will be espe- 
cially problematic. For example, GRACE will not sample a complete phase cycle 
of Ki until nearly four years have passed, and such long periods will contain 
substantial power from non-tidal signals. The aliasing of certain tides into semi- 
annual and near-annual periods will also make solar tide improvements difficult. 
Lunar tides will be less problematic, and improvements even to M 2 will be most 
welcome. 

Space limitations have forced us to concentrate on the subject of tides as a 
correction to new satellite gravimetry missions. But if history is a guide, then 
the markedly enhanced precisions of these new measurements will also yield new 
information about the tides — solid and atmospheric as well as oceanic — and such 
improvements can in turn yield further information about other aspects of the Earth 
system. As just one example, improvements in degree-2 and degree-3 tides shed 
light on the Earth’s anelasticity at frequencies well outside the seismic band. De- 
termination of these low-degree coefficients should improve significantly with the 
new satellite data. 
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